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GENERALIZED TWO-STAGE DATA ESTIMATION 
[00011 FIELD OF INVENTION 

[0002] The present invention relates to wireless communication systems. 
More particularly, the present invention is directed to data estimation in such 
systems. 

[0003] BACKGROUND 

[0004] In wireless systems, joint detection (JD) is used to mitigate inter- 
symbol interference (ISD and multiple-access interference (MAI). JD is 
characterized by good performance but high complexity. Even using approximate 
Cholesky or block Fourier transforms with Cholesky decomposition algorithms, 
the complexity of JD is still very high. When JD is adopted in a wireless receiver, 
its complexity prevents the receiver from being implemented efficiently. This 
evidences the need for alternative algorithms that are not only simple in 
implementation but also good in performance. 

[0005] To overcome this problem, prior art receivers based on a channel 
equalizer followed by a code despreader have been developed. These types of 
receivers are called single user detection (SUD) receivers because, contrary to JD 
receivers, the detection process does not require the knowledge of channelization 
codes of other users. SUD tends to not exhibit the same performance as JD for 
most data rates of interest, even though its complexity is very low. Accordingly, 
there exists a need for low complexity high performance data detectors. 

[0006] SUMMARY 

[0007] Symbols are to be recovered from signals received in a shared 
spectrum. Codes of the signals received in the shared spectrum are processed 
using a block Fourier transform (FT), producing a code block diagonal matrix. A 
channel response of the received signals is estimated. The channel response is 
extended and modified to produce a block circulant matrix and a block FT is 
taken, producing a channel response block diagonal matrix. The code block 
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diagonal matrix is combined with the channel response block diagonal matrix. 
The received signals are sampled and processed using the combined code block 
diagonal matrix and the channel response block diagonal matrix with a Cholesky 
algorithm. A block inverse FT is performed on a result of the Cholesky algorithm 
to produce spread symbols. The spread symbols Eire despread to recover symbols 
of the received signals. 

[0008] BRIEF DESCRIPTION OF THE DRAWINGS 

[0009] Figure 1 is a block diagram showing a two stage data detection. 

[0010] Figure 2 is a block diagram of an embodiment of two-stage data 

detection. 

[0011] Figure 3 is a block diagram of code assignment to reduce the 
complexity of two-stage data detection. 

[0012] Figure 4A-4D are block diagrams of utilizing look-up tables to 
determine ^R. 



[0013] DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT(S) 
[0014] The present invention will be described with reference to the 
drawing figures where like numerals represent like elements throughout. 
[0015] A two stage data estimator can be used in WTRU or base station, 
when all of the communications to be detected by the estimator experience a 
similar channel response. Although the following is described in conjunction 
with the preferred proposed third generation partnership project (3GPP) 
wideband code division multiple access (W-CDMA) communication system, it is 
applicable to other systems. 

[0016] Figure 1 is a simplified block diagram of a receiver using a two stage 
data estimator 55. An antenna 50 or antenna array receives radio frequency 
signals. The signals are sampled by a sampling device 51, typically at the chip 
rate or at a multiple of the chip rate, producing a received vector r . A channel 
estimation device 53 using a reference signal, such as a midamble sequence or 
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pilot code, estimates the channel response for the received signals as a channel 
response matrix H. The channel estimation device 53 also estimates the noise 
variance, o 2 . 

[0017] The channel equalizer 52 takes the received vector r and equalizes 
it using the channel response matrix H and a the noise variance o 2 , producing a 
spread symbol vector s . Using codes C of the received signals, a despreader 54 
despreads the spread symbol vectors , producing the estimated symbols d . 
[0018] With joint detection (JD), a minimum mean square error (MMSE) 
formula with respect to the symbol vector - can be expressed as: 

Equation (1) 

or 

d = R d A H {AR d A H +R n T x Li 

Equation (2) 

£ is the estimate of ^ , £ is the received signal vector, A is the system matrix, R » 

is the covariance matrix of noise sequence, is the covariance matrix of the 

symbol sequence and the notation denotes the comply conjugate transform 
(Hermitian) operation. The dimensions and structures of the above vectors and 
matrixes depend on specific system design. Usually, different systems have 
different system parameters such as frame structure, length of data field and 
length of delay spread. 

[0019] The matrix A has the different values of dimensions for different 
systems and the dimensions of matrix A depend on the length of data field, 
number of codes, spreading factor and length of delay spread. By way of 
example, for the transmission of 8 codes with spreading factor 16 each, the 
matrix A has dimensions of 1032 by 488 for a WCDMA TDD system if burst type 
1 is used and for delay spread of 57 chips long, while matrix A has dimensions of 
367 by 176 for TD-SCDMA system for a delay spread of 16 chips long. 
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[0020] Assuming white noise and uncorrected symbols with unity energy, 

K - a% I and = ' , where I denotes the identity matrix. Substitution of these 
into Equations 1 and 2 results in: 
d={A H A + o 2 iy x A H r 7 

Equation (3) 

or 

d = A H (AA H +a 2 iy l r. 

Equation (4) 

[0021] The received signal can be viewed as a composite signal, denoted by 
1 9 passed through a single channel. The received signal £ may be represented by 
I - Hi, where H is the channel response matrix and £ is the composite spread 
signal. H takes the form of: 

ft, fto 
. ft, . 

. ft 0 
ft, 

Equation (5) 

[0022] In Equation (5), W is the length of the channel response, and is 
therefore equal to the length of the delay spread. Typically W=57 for W-CDMA 
time division duplex (TDD) burst type 1 and W=16 for time division synchronous 
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CDMA (TD-SCDMA). The composite spread signal £ can be 
where the symbol vector ^ is: 

and the code matrix C is: 

c=[c"\c«\...,c {K) ] 



with: 



Equation (6) 



Equation (7) 



Equation (8) 



[0023] Q, K and N * denote the spread factor (SF), the number of active 
codes and the number of symbols carried on each channelization code, 

respectively. c '* } is the i* element of the ** code. The matrix C is a matrix of 
size N.-Q^N.-K. 

[0024] Substitution of A = HC into Equation (4) results in: 
d = C H H H (HR c H H +a 2 ir l r 

Equation (9) 



[0025] 
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[0026] Rc = CC ♦ If - denotes the estimated spread signal, Equation (9) 
can be expressed in two stages: 
[0027] Stage 1: 

s = H H (HR c H H +v 2 iy l r 

Equation (10) 

[0028] Stage 2: ' 
d = C H S 

Equation (11) 

[0029] The first stage is the stage of generalized channel equalization. It 

estimates the spread signal £ by an equalization process per Equation 10. The 

second stage is the despreading stage. The symbol sequence ^ is recovered by 
a despreading process per Equation 11. 

[0030] The matrix R * in Equation 9 is a block diagonal matrix of the 
form: 

*o 

O 



Equation (12) 

[0031] The block R o in the .diagonal is a square matrix of size Q. The 
matrix R c is a square matrix of size ^* ' ^ . 

[0032] Because the matrix R c is a block circular matrix, the block Fast Fourier 

R 

transform (FFT) can be used to realize the algorithm. With this approach the matrix c can 
be decomposed as: 

R c = F tQ) A * F iQ) 

Equation (13) 

with 
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Equation (14) 

[0033] f n* is the N * -point PFT matrix, 7 Q is the identity matrix of size Q 
and the notation ® is the Kronecker product. By definition, the Kronecker 
product Z of matrix X and Y, (Z = X ® 7 ) is: 

'^y *, 2 r a x xn y~ 



z = 



M O 
x My Y x MX Y x MN Y 



Equation (15) 

x m is the ( m > n ) th element of matrix X. For each ^(fl) , a Ns-pointEFT is performed Q times. 
A/? is a block-diagonal matrix whose diagonal blocks are: 
F (C) /e c 0a:e) That is, 

diag(A R ) = F iQ) R c (:*:Q) 9 

Equation (16) 

^cC' 1 : 2) denotes the first Q columns of matrix^ . 

[0034] The block circular matrix can be decomposed into simple and 
efficient FFT components, making a matrix inverse more efficient and less 
complex. Usually, the large matrix inverse is more efficient when it is 
performed in the frequency domain rather than in a time domain. For this 
reason, it is advantage to use FFT and the use of a block circular matrix 
enables efficient FFT implementation. With proper partition, the matrix H 
can be expressed as a approximate block circular matrix of the form: 
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H. H n 



H 



L-l 



L-l 



*0 

H 



Equation (17) 

where each ff ' , 1 = 0>1 " L-1 is a square matrix of size Q. L is the number of 
data symbols affected by the delay spread of propagation channel is expressed 



as: 



Q 



[0035] ^ - Equation (18) 

[0036] To enable block FFT decomposition, H can be extended and modified 
into an exactly block circular matrix of the form: 



[0037] 



H 



L-l 



Bo 



H 



L-l 



H 



L-l 



H 2 H x H 0 



Equation (19) 

[0038] The block circular matrix H c is obtained by expanding the columns 
of matrix H in Equation (17) by circularly down-shifting one element block 
successively. 
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[0039] The matrix *^can be decomposed by block PFT as: 

H C =zF <Q) A H F iQ) 

Equation (20) 

[0040] A n is a block-diagonal matrix whose diagonal blocks are 

F m H c M : Q) as dwg(A ff ) = F {Q) H C (:1 : 0 

Equation (21) 

[0041] H C (: 9 1:Q) denotes the first Q columns of matrix H c. From 

Equation (20), H" can be defined as 

H c = F iQ) AH H F iQ) Equation (22) 

Substituting matrix R c and H c into Equation 10, £ is obtained: 
[0042] l=^(A H AX +cr 2 ir l F {Q) r 

Equation (23) 

[0043] For a zero forcing (ZF) solution, equation 19 is simplified to 

1= F (Q) A ~R^~H F {Q)L 

Equation (24) 

[0044] The matrix inverse in Equations (23) and (24) can be performed 
using Cholesky decomposition and forward and backward substitutions. 
[0045] In a special case of K=SF, where (the number of active codes 

equals the spreading factor), the matrix R c becomes a scalar-diagonal matrix 
with identical diagonal elements equal to SF. In this case, Equations (10) and 
(11) reduce to: 

s = H H (HH H +^iy x L 

Equation (25) 

and 

Equation (26) 

[0046] Equation (25) can also be expressed in the form of: 
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$ = (H H H+£iy l H H r 

Equation (27) 

[0047] With FFT, Equations (25) and (27) can be realized by: 
S = F- l A* H (A H A* H +$ir l Fr 

Equation (28) 

and 

Equation (29) 

respectively. A " is a diagonal matrix whose diagonal is F ' H ( : X) in which 
#M) denotes the first column of matrix H. The notation (.)* denotes conjugate 
operator. 

[0048] Figure 2 is a preferred block diagram of the channel equalizer 15. 
A code matrix C is input into the channel equalizer 15. A Hermitian device 30 
takes a complex conjugate transpose of the code matrix C, C H . The code matrix 
C and its Hermitian are multiplied by a multiplier 32, producing CC H . A block 
FT performed on CC H ' producing block diagonal matrix A* . 
[0049] The channel response matrix H is extended and modified by an 
extend and modify device 36, producing fl°. A block FT 38 takes B° and 
produces block diagonal matrix A H . A multiplier multiplies A H and A* 
together, producing A H A R . A Hermitian device 42 takes the complex conjugate 
transpose of A H , producing A^. A multiplier 44 multiplies A^ to A H A Rf 
producing A H A R A* to o 2 I , A producing A H A R A H H + cr 2 / . 
[0050] A Cholesky decomposition device 48 produces a Cholesky factor. A 
block FT 50 takes a block FT of the received vector r. Using the Cholesky 
factor and the FT of r , forward and backward substitution are performed by a 
forward substitution device 52 and backward substitution device 54. 
[005 1] A conjugation device 56 takes the conjugate of A H , producing A „ . 
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The result of backward substitution is multiplies to A* H . A block inverse FT 60 

A 

takes a block inverse FT of the multiplied result, producing s. 
[0052] According to another embodiment of the present invention, an 
approximate solution is provided in which the generalized two-stage data 
detection process is a block-diagonal-approximation. The block-diagonal- 
approximation includes off-diagonal entries as well as the diagonal entries in 
the approximation process. 

[0053] As an example, the case of four channelization codes is considered. 
Ro. a combination of four channelization codes, comprises a constant block 
diagonal part, which does not vary with the different combinations of the codes, 
and an edge part which changes with the combinations. In general Ro has the 
structure of: 



[0054] 
[0055] 



c 
c 

X 
X 



X X 
X X 



c 
c 

X 
X 



Equation (30) 

where elements denoted as c represent constants and are always 
equal to the number of channelization codes, i.e., c = K m The elements 
designated as x represent some variables whose values and locations vary with 
different combinations of channelization codes. Their locations vary following 
certain patterns depending on combinations of codes. As a result only a few of 
them are non-zero. When code power is considered and is not unity power, the 
element c equals the total power of transmitted codes. A good approximation of 
the matrix Rois to include the constant part and ignore the variable part as: 



-11- 



WO 2004/064298 



PCT/US2004/000463 



[0056] 



c 
c 



0 



c c 
c c 



Equation (31) 

[0057] In this case, the approximation ^> contains only a constant part. *° 
depends only on the number of active codes regardless of which codes are 

transmitted, and *c can be decomposed as shown is Equation (13). The block 

diagonal of A R or F (fiA( : > 1: 2) can be pre-calculated using an FFT for different 
numbers of codes and stored as a look-up table. This reduces the computational 

complexity by not computing F (0) R c( : ^'-Q\ In the case, that code power is 
considered and is not unity power, the element c becomes total power of active 

codes, (i.e., c = F r in which p t is the total power of active codes). The matrix *o 
can be expressed as 



^0 — P avs 



K K 
K K 



K K 
K K 



K K 
K K 



[0058] 



K K 
K K 



Equation (32) 
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[0059] where is the average code power obtained by ^ K . In this 
case, a scaling Pm * should be applied in the process. 

[0060] Other variants of block-diagonal approximation method can be 
derived by including more entries other than the constant block-diagonal part. 
This improves performance but entails more complexity because by including 

variable entries the FFT for F iQ) R c^ 1 : 2) has to be now recalculated as needed 
if the codes change. The use of more entries enhances the exact solution as all 
of the off-diagonal entries are included for processing. 

[006 1] At a given number of channelization codes, one can derive the code 
sets for different combinations of channelization codes that have common 
constant part of the correlation matrix whose values are equal to the number of 
channelization codes, or the total power of channelization codes when the code 
does not have unity code power. To facilitate the low complexity 
implementation, the assignment of channelization codes or resource units can 
be made following the rules that a code set is randomly picked among the code 
sets that have common constant part and those codes in the picked code set are 
assigned. For example of assignment of four codes, the code sets [1,2,3,4], 
[5,6,7,8], [9,10,11,12], .... have the common constant part in their correlation 
matrix. When channel assignment of four codes is made, one of those code sets 
should be used for optimal computational efficiency. 

[0062] Figure 3 is a flow diagram of such a channel code assignment. 
Codes sets having a constant part are determined, step 100. When assigning 
codes, the code sets having the constant part are used, step 102. 
[0063] Figures 4A, 4B, 4C and 4D are illustrations of preferred circuits for 
reducing the complexity in calculating A R . In Figure 4A, the number of codes 
processed by the two stage data detector are put in a look-up table 62 and the 
A R associated with that code number is used. In Figure 4B, the number of 
codes processed by the two stage data detector are put in a look-up table 62 and 
an unsealed A R is produced. The unsealed A R is scaled, such as by a multiplier 
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66 by P BVB , producing A R . 

[0064] In Figure 4C, the code matrix C or code identifier is input into a 
look-up table 68. Using the look-up table 68, the A* is determined. The 
unsealed A R is scaled, such as by a multiplier 72 by P avg , producing A* . In 
Figure 4D, the code matrix C or code identifier is input into a look-up table 70, 

producing an unsealed A*. 

* * * 
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CLAIMS 

What is claimed is: 

1. A method for recovering symbols from signals received in a shared 
spectrum, the method comprising: 

processing codes of the signals received in the shared spectrum using a 
block Fourier transform (FT) and producing a code block diagonal matrix; 

estimating a channel response of the received signals; 

extending and modifying the channel response to produce a block circulant 
matrix and taking a block FT and producing a channel response block diagonal 
matrix; 

combining the code block diagonal matrix and the channel response block 
diagonal matrix; 

sampling the received signals; 

processing the received signal using the combined code block diagonal 
matrix and the channel response block diagonal matrix with a Cholesky 
algorithm; 

performing a block inverse FT on a result of the Cholesky algorithm to 
produce spread symbols; and 

despreading the spread symbols to recover symbols of the received signals. 

2. The method of claim 1 wherein the Cholesky algorithm includes 
determining a Cholesky factor and performing forward and backward 
substitution. 
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3. The method of claim 1 wherein the combining the code block 
diagonal matrix and the channel response block diagonal matrix includes adding 
factor of the noise variance multiplied with an identity matrix. 

4. The method of claim 1 wherein the code block diagonal matrix is 
produced by multiplying a code matrix with a complex conjugate transpose of the 
code matrix and taking a block FT of a result of the multiplying. 

5. The method of claim 1 wherein the code block diagonal matrix is 

j 

produced by inputting a number of codes of interest into a look-up table. 

6. The method of claim 1 wherein the code block diagonal matrix is 
produced by inputting a number of codes of interest into a look-up table and 
scaling a resulting diagonal block matrix from the look-up table by an average 
power level. 

7. The method of claim 1 wherein the code block diagonal matrix is 
produced by inputting code identifiers of the received signals into a look-up table. 

8. The method of claim 1 wherein the code block diagonal matrix is 
produced by inputting code identifiers of the received signals into a look-up table 
and scaling a resulting diagonal block matrix from the look-up table by an 
average power level. 
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9. The method of claim 1 wherein the code block diagonal matrix is 
produced by inputting codes of the received signal into a look-up table. 



10. The method of claim 1 wherein the code block diagonal matrix is 
produced by inputting inputing codes of the received signal into a look-up table 
and scaling a resulting diagonal block matrix from the look-up table by an 
average power level. 

11. A wireless transmit/receive unit (WTRU) for use in recovering 
symbols from signals received in a shared spectrum, the WTRU comprising: 

means for processing codes of the signals received in the shared spectrum 
using a block Fourier transform (FT) and producing a code block diagonal matrix; 

means for estimating a channel response of the received signals; 

means for extending and modifying the channel response to produce a 
block circulant matrix and taking a block FT and producing a channel response 
block diagonal matrix; 

means for combining the code block diagonal matrix and the channel 1 
response block diagonal matrix; 

means for sampling the received signals; 

means for processing the received signal using the combined code block 
diagonal matrix and the channel response block diagonal matrix with a Cholesky 
algorithm; 

means for performing a block inverse FT on a result of the Cholesky 
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algorithm to produce spread symbols; and 

means for despreading the spread symbols to recover symbols of the 
received signals. 

12. The WTRU of claim 11 wherein the Cholesky algorithm includes 
determining a Cholesky factor and performing forward and backward 
substitution. 

13. The WTRU of claim 11 wherein the combining the code block 
diagonal matrix and the channel response block diagonal matrix includes adding 
factor of the noise variance multiplied with an identity matrix. 

14. The WTRU of claim 11 wherein the code block diagonal matrix is 
produced by multiplying a code matrix with a complex conjugate transpose of the 
code matrix and taking a block FT of a result of the multiplying. 

15. The WTRU of claim 11 wherein the code block diagonal matrix is 
produced by inputting a number of codes of interest into a look-up table. 

16. The WTRU of claim 11 wherein the code block diagonal matrix is 
produced by inputting a number of codes of interest into a look-up table and 
scaling a resulting diagonal block matrix from the look-up table by an average 
power level. 
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17. The WTRU of claim 11 wherein the code block diagonal matrix is 
produced by inputting code identifiers of the received signals into a look-up table. 

18. The WTRU of claim 11 wherein the code block diagonal matrix is 
produced by inputting code identifiers of the received signals into a look-up table 
and scaling a resulting diagonal block matrix from the look-up table by an 
average power level. 

19. The WTRU of claim 11 wherein the code block diagonal matrix is 
produced by inputting codes of the received signal into a look-up table. 

20. The WTRU of claim 11 wherein the code block diagonal matrix is 
produced by inputting inputing codes of the received signal into a look-up table 
and scaling a resulting diagonal block matrix from the look-up table by an 
average power level. 

31. A wireless transmit/receive unit (WTRU) for use in recovering 
symbols from signals received in a shared spectrum, the WTRU comprising: 

a block Fourier transform (FT) device for processing codes of the signals 
received in the shared spectrum using a block FT and producing a code block 
diagonal matrix; 

a channel estimation device for estimating a channel response of the 
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received signals; 

an extending and modifying block for extending and modifying the channel 
response to produce a block circulant matrix and taking a block FT and producing 
a channel response block diagonal matrix; 

a circuit for combining the code block diagonal matrix and the channel 
response block diagonal matrix; 

a sampling device for sampling the received signals; 

a Cholesky decomposition device and forward and backward substitution 
devices for processing the received signal using the combined code block diagonal 
matrix and the channel response block diagonal matrix with a Cholesky 
algorithm; 

an inverse block FT device for performing a block inverse FT on an output 
of the backward substitution device to produce spread symbols; and 

a despreader for despreading the spread symbols to recover symbols of the 
received signals, 

32. The WTRU of claim 31 wherein the circuit for combining comprises 
a two multipliers. 

33. The WTRU of claim 31 wherein the circuit for co m bining includes 
adding factor of the noise variance multiplied with an identity matrix. 

34. The WTRU of claim 3 1 further comprising a Hermitian device and a 
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multiplier for multiplying a code matrix with a complex conjugate transpose of 
the code matrix. 



35. The WTRU of claim 31 further comprising a look-up table wherein 
the code block diagonal matrix is produced by inputting a number of codes of 
interest into a look-up table. 

36. The WTRU of claim 31 comprising a look-up table and a multiplier 
wherein the code block diagonal matrix is produced by inputting a number of 
codes of interest into a look-up table and multiplying a resulting diagonal block 
matrix from the look-up table by an average power level. 

37 . The WTRU of claim 3 1 further comprising a look-up table wherein 
the code block diagonal matrix is produced by inputting code identifiers of the 
received signals into a look-up table. 

38. The WTRU of claim 31 comprising a look-up table and a multiplier 
wherein the code block diagonal matrix is produced by inputting code identifiers 
of the received signals into a look-up table and multiplying a resulting diagonal 
block matrix from the look-up table by an average power level. 

39. The WTRU of claim 31 further comprising a look-up table wherein 
the code block diagonal matrix is produced by inputting codes of the received 
signal into a look-up table. 
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40. The WTRU of claim 31 comprising a look-up table and a multiplier 
wherein the code block diagonal matrix is produced by inputting inputting codes 
of the received signal into a look-up table and multiplying a resulting diagonal 
block matrix from the look-up table by an average power level. 

51. A base station for use in recovering symbols from signals received in 
a shared spectrum, the base station comprising: 

means for processing codes of the signals received in the shared spectrum 
using a block Fourier transform (FT) and producing a code block diagonal matrix; 

means for estimating a channel response of the received signals; 

means for extending and modifying the channel response to produce a 
block circulant matrix and taking a block FT and producing a channel response 
block diagonal matrix; 

means for combining the code block diagonal matrix and the channel 
response block diagonal matrix; 

means for sampling the received signals; 

means for processing the received signal using the combined code block 
diagonal matrix and the channel response block diagonal matrix with a Cholesky 
algorithm; 

means for performing a block inverse FT on a result of the Cholesky 
algorithm to produce spread symbols; and 

means for despreading the spread symbols to recover symbols of the 
received signals. 
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52. The base station of claim 51 wherein the Cholesky algorithm 
includes determining a Cholesky factor and performing forward and backward 
substitution. 

53. The base station of claim 51 wherein the combining the code block 
diagonal matrix and the channel response block diagonal matrix includes adding 
factor of the noise variance multiplied with an identity matrix. 

54. The base station of claim 51 wherein the code block diagonal matrix 
is produced by multiplying a code matrix with a complex conjugate transpose of 
the code matrix and taking a block FT of a result of the multiplying. 

55 . The base station of claim 5 1 wherein the code block diagonal matrix 
is produced by inputting a number of codes of interest into a look-up table. 

56. The base station of claim 51 wherein the code block diagonal matrix 
is produced by inputting a number of codes of interest into a look-up table and 
scaling a resulting diagonal block matrix from the look-up table by an average 
power level. 

57. The base station of claim 51 wherein the code block diagonal matrix 
is produced by inputting code identifiers of the received signals into a look-up 
table. 
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58 . The base station of claim 5 1 wherein the code block diagonal matrix 
is produced by inputting code identifiers of the received signals into a look-up 
table and scaling a resulting diagonal block matrix from the look-up table by an 
average power level. 

59. The base station of claim 5 1 wherein the code block diagonal matrix 
is produced by inputting codes of the received signal into a look-up table. 

60. The base station of claim 51 wherein the code block diagonal matrix 
is produced by inputting inputting codes of the received signal into a look-up 
table and scaling a resulting diagonal block matrix from the look-up table by an 
average power level. 

61. A base station for use in recovering symbols from signals received in 
a shared spectrum, the base station comprising: 

a block Fourier transform (FT) device for processing codes of the signals 
received in the shared spectrum using a block FT and producing a code block 
diagonal matrix; 

a channel estimation device for estimating a channel response of the 
received signals; 

an extending and modifying block for extending and modifying the channel 
response to produce a block circulant matrix and taking a block FT and producing 
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a channel response block diagonal matrix; 

a circuit for combining the code block diagonal matrix and the channel 
response block diagonal matrix; 

a sampling device for sampling the received signals; 

a Cholesky decomposition device and forward and backward substitution 
devices for processing the received signal using the combined code block diagonal 
matrix and the channel response block diagonal matrix with a Cholesky 
algorithm; 

an inverse block FT device for performing a block inverse FT on an output 
of the backward substitution device to produce spread symbols; and 

a despreader for despreading the spread symbols to recover symbols of the 
received signals. 

62. The base station of claim 61 wherein the circuit for combining 
comprises a two multipliers. 

63. The base station of claim 61 wherein the circuit for combining 
includes adding factor of the noise variance multiplied with an identity matrix. 

64. The base station of claim 61 further comprising a Hermitian device 
and a multiplier for multiplying a code matrix with a complex conjugate 
transpose of the code matrix. 
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65. The base station of claim 61 further comprising a look-up table 
wherein the code block diagonal matrix is produced by inputting a number of 
codes of interest into a look-up table. 

66. The base station of claim 61 comprising a look-up table and a 
multiplier wherein the code block diagonal matrix is produced by inputting a 
number of codes of interest into a look-up table and multiplying a resulting 
diagonal block matrix from the look-up table by an average power level. 

67. The base station of claim 61 further comprising a look-up table 
wherein the code block diagonal matrix is produced by inputting code identifiers 
of the received signals into a look-up table. 

68. The base station of claim 61 comprising a look-up table and a 
multiplier wherein the code block diagonal matrix is produced by inputting code 
identifiers of the received signals into a look-up table and multiplying a resulting 
diagonal block matrix from the look-up table by an average power level. 

69. The base station of claim 61 further comprising a look-up table 
wherein the code block diagonal matrix is produced by inputting codes of the 
received signal into a look-up table. 

70. The base station of claim 61 comprising a look-up table and a 
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multiplier wherein the code block diagonal matrix is produced by inputting 
inputting codes of the received signal into a look-up table and multiplying a 
resulting diagonal block matrix from the look-up table by an average power level. 
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